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Abstract 

The infinite limit of Matrix Theory in 4 and 10 dimensions is described in terms of Moyal 
Brackets. In those dimensions there exists a BogomoPnyi bound to the Euclideanized version 
of these equations, which guarantees that solutions of the first order equations also solve the 
second order Matrix Theory equations. A general construction of such solutions in terms 
of a representation of the target space co-ordinates as non-local spinor bilinears, which are 
generalisations of the standard Wigner functions on phase space, is given. 



1 Introduction 

The purpose of the present paper is to clarify the meaning of the large N limit in M-theory and 
discuss the features of the formalism peculiar to the situation of 2 and 8 transverse dimensions. 
When this limit is approached by way of the Moyal Bracket formalism, the matrix M-theory takes on 
an aspect resembling a generalisation of the Moyal formulation of Quantum Mechanics ||14|| , in terms 
of Wigner phase space distributions. The target space co-ordinates in 4 and 10-dimensions will be 
represented in terms of matrix elements between spinor wave functions as Wigner distributions. 
0] The non-local character of M-theory [[| is perhaps reflected in this construction. The paper 
begins with a short review of the usual interpretation of the large N limit before developing the 
interpretation in terms of Moyal Brackets. The essential simplicity of the theory in 4 and 10 
dimensions resides in the fact that the Euclidean form of the Lagrangian can be written as a sum 
of squares of first order terms with positive relative signs, whose vanishing implies the second order 
equations of motion, thus giving rise to BPS states. These first order equations may be interpreted 
as a generalisation of the Liouville equation for distribution functions. This is the basis of the 
suggested analogy with Quantum Mechanics. This situation is analysed in detail to give a general 
principle for the construction of solutions to these equations. The 10- dimensional situation is 
complicated by the presence of equations of constraint. In the 4-dimensional case these constraints 
are absent, and an alternative solution procedure is outlined. 

1 On research leave from the University of Durham, U.K. 
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2 Setting the Scene 

The Lagrangian describing matrix models has been described some time ago ||. In recent months 
the paper of Banks et al for a matrix model of M-theory || has stimulated the production of a 
profusion of papers elaborating the issue 0. We take one such paper, as a convenient reference 
point. The Lagrangian to be used, taken directly from takes the form of a two-dimensional 
M = 8 supersymmetric U(N) Yang-Mills theory with the action 

S = ^7 / Tr ((ZVO 2 + T p + g\Fl v - ±[X\ X^ + jfli[X\ *]) dadr. (1) 

Here the 8 scalar fields X % are N x N hermitian matrices, as are the 8 fermionic fields and 8^. 
The fields X 1 , 8 a , 8 a transform respectively as the 8^ vector, and 8 S and 8 C spinor representations of 
the SO (8) R-symmetry group of transversal rotations. The two dimesional world sheet is a cylinder 
parametrised by co-ordinates a, r, with < a < 2n. This article goes on to identify vacuum 
configurations of the theory in the infra-red limit (g s \— > 0) with matrices unitarily equivalent to 
diagonal ones i.e. 

X» = U~ l x^U (2) 

where x^ are diagonal matrices with entries x^. In the simplest description of the model each 
of the eigenvalues Xu, 8%, 8% is considered to represent a Green-Schwarz string in the light cone 
gauge, and the entire model is described in this limit as a gas of N such strings. In a more general 
situation envisaged in (|T]) as a, the parameter in the compactified direction increases from to 2tt, 
instead of matching xu(0) with xh(2ti) several subsets of the eigenvalues may become permuted 
and the eigenvalues return to their original values (for a given subset) only after m k turns round 
the compact direction. In this case the theory describes yj N k strings, where 

J2m k N k = N. (3) 

The authors of ([I]) then proceed to argue that the bulk of the energy in the infinite momentum 
frame will be carried by long strings in the limit of very large N. There is, however, another way to 
consider the infinite limit which ends with a theory formulated in terms of Poisson Brackets. This 
theory possesses some unusual features; the bosonic sector possesses hidden invariances due to Bars 
[§] which apparently are lost in the finite matrix version, and which represent the remnants of volume 
preserving diffeomorphisms. Furthermore, when the transverse space is either 2 or 8-dimensional 
the first order (bosonic) equations of motion can be written which provide a Bogomol'nyi bound on 
the action. These first order equations furthermore, also possess the additional symmetry. It is to 
be expected that something similar happens in the case of 4 transverse dimensions, but this case is 
not analysed in the present article. Another approach to BPS configurations can be found in |1. 



3 Passage to the Infinite Limit 

The passage to the iV — > oo limit can be conceived as taking place in two steps; First of all the fields 
are described by functions of two 'phase space' variables a, (3 as well as cr, r which parametrise the 
world sheet instead of matrices, and expressions of the form J Tr[X*, 2 da by 

J (- sm{X\ X ] }) 2 dad(3do (4) 
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where sin{X*, X J } is the sine, or Moyal Bracket, ||, with deformation parameter A which is defined 
as the imaginary part of the star product (*); 

X»*X U = lim e iA( ^- a ^ } X^(a, P, cr)X u (a\ (3', a). (5) 

a 1 — >a 



2tt 



The point of this construction is that in the limiting points A — > — the Moyal brackets [10 



re- 



produce the commutators of N x N matrices through the association of the components X^ n of a 
matrix X^ with the Fourier modes of a function X^{a, (3, a), periodic in a, (5. (Strictly speaking, 
N should be odd for this association to work.) The fermionic terms in ([!]) are replaced by the cosine 
bracket, i.e. the real part of the star product. Thus the action becomes 

S M b = 7^7 J {{D^Xf + cos{# T , p} + g 2 s TrF 2 u ) dadpdadr 

~ I \(-T^sm{X\X j }) 2 - — cosi^, jsinpf\0}} ) dadpdadr. (6) 
J \ ^9; 9s A J 

The final step involves letting N grow indefinitely, or equivalently, letting A — > 0. The final form of 
the action S is then expressed in terms of ordinary Poisson Brackets; 

Spb = ^7 / {{D,X f + 9 T , p + g 2 TrF^ u - (±{X\ X^ dad(3dadr 

+ J (^- dT Hi 0}\ dadpdadr. (7) 

In this interpretation the theory describes a 3-membrane wrapped around the compactified direction. 
In order to treat the system dynamically, the longitudinal and timelike modes should be re-instated. 



4 Membranes in 8 transverse dimensions 

A remarkable feature of the situation with 8 transverse directions is that the theory admits a 
class of solutions obtainable from a first order formulation, thanks to the existence of a a self- 
dual (antisymmetric) 4-tensor T^ vpa in 10-dimensions which is an extension of that discovered in 
8-dimensions [JO]] and which is an analogue of the 4-dimensional fully antisymmetric tensor e pvpa . 



The equations quoted below are of similar form to those exhibited as equations (63) in [fL2jl , but 
with a more physically appropriate reinterpretation of the nature of the variables. X° describes 
the timelike directions, X 9 the longitudinal and X 1 , i = 1 . . . 8 the transverse directions. We shall 
also ignore the dependence upon r, which is usually interpreted as world-sheet time and treat the 
additional independent variable as a . This procedure is tanatamount to working in a Euclidean, 
rather than Lorentzian space, as it is characteristic of equations of self-duality that they should be 
formulated in a Euclidean space to ensure that solutions will be real. This is the same sign as that 



adopted in [13 



Then the 10-dimensional membrane has a Lagrangian density in the form of an SU(oo) pure 
Yang Mills theory (0). The situation is further simplified to one of dependence upon only one of 
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the variables a (apart from the a, (3 dependence of the gauge potentials X^(t, a, (3)). We work in 
a gauge where X° = constant. The Lagrangian density is 



C 



\{d a x»Y + \{X^ x»y 



where curly brackets denote the Poisson Bracket with respect to the variables a, /3, i.e. 

dX» dX v dX» dX v 



{X", X u } 



da 8(3 8(3 da 



(8) 



(9) 



and the theory is equivalent to an SU(oo) Yang Mills with dependence upon only one variable. It 
is perhaps worth recalling that in the case of finite X Savvidy and others have shown that the 
classical second order equations for such a theory in the case when that variable is time exhibit 
chaotic solutions It is not clear whether this statement persists when the large N limit is 

taken. 

The first order set of equations is then given by 







d^X 1 + {X 2 , 


X 9 } 


= 






d a x 2 + {x 9 , 


x 1 } 


= 






d a X 3 + {X 4 , 


X 9 } 


= 






d a X 4 + {X 9 , 


X 3 } 


= 






d a X 5 + {X 6 , 


X 9 } 


= 






d a X 6 + {X 9 , 


X 5 } 


= 






d a X 7 + {X 8 , 


X 9 } 


= 






<9 CT X 8 + {X 9 , 


X 7 } 


= 


d a X 9 + {X 1 , 


X 2 } + {X 3 , 


X 4 } + {X 5 , X 6 } + {X 7 , 


X 8 } 


= 


{x\ 


X 3 } + {X 4 , 


X 2 } + {X 5 , X 7 } + {X 8 , 


X 6 } 


= 


{x\ 


X 4 } + {X 2 , 


X 3 } + {X s , X 5 } + {X 7 , 


X 6 } 


= 


{x\ 


X 5 } + {X 4 , 


X 8 } + {X 7 , X 3 } + {X 6 , 


X 2 } 


= 


{x\ 


X 6 } + {X 2 , 


X 5 } + {X 3 , X 8 } + {X 4 , 


X 7 } 


= 


{x\ 


X 7 } + {X 3 , 


X 5 } + {X s , X 2 } + {X 6 , 


X 4 } 


= 


{x\ 


X 8 } + {X 5 , 


X 4 } + {X 2 , X 7 } + {X 6 , 


X 3 } 


= 



(10) 

It is easy to verify that the second order equations coming from the Lagrangian @ are satisfied by 
solutions of this system and that the sums of squares of these equations give the Lagrangian density 
up to divergences whether the brackets {} are taken to signify commutators, Moyal Brackets or 
Poisson brackets. Thus the solutions provide a bound on the action. In the case of Poisson brackets, 
in fact only the cross terms involving a a derivative give divergences; the others vanish identically 
in virtue of what we call the 'Saturday Identity', 



{/, g}{h, k} + {/, h}{k, g} + {/, k}{g, h} = 0, 



which holds for Poisson Brackets on a 2-dimensional phase-space-but not for commutators, nor for 
the Moyal Bracket form. These equations will give ries to solitonic solutions. 



5 



What we have in ( |10D in the case where X° = constant is a set of coupled Nahm equations; 
When expressed in terms of complex combinations defined by 

Zl = X 1 + iX 2 - z 2 =X 3 + iX*; z 3 = X 5 + iX 6 ; Zl = X 7 + iX s ; H = X 9 , (11) 

these take the following form 



d a H + {zi, z x } + {z 2 , 



d a z\ + i{H, zi} 


= 


d a z 2 + i{H, z 2 } 


= 


d a z 3 + i{H, z 3 } 


= 


d a z 4 + i{H, z 4 } 


= 


z 2 } + {23, z 3 } + {z 4 , z 4 } 


= 


{zi, z 2 } + {z 3 , z 4 } 


= 


{z x , z 3 } + {z 4 , z 2 } 


= 


{zi, z 4 } + {z 2 , z 3 } 


= 



5 Steps towards a solution 

The main purpose of writing the equations in complex form is to emphasise the similarity of the 
equations to the phase space formulation of Quantum Mechanics. 

The observation that the first four equations are of the (Quantum) Liouville type provides a 
clue towards a possible form of solution. To emphasise the similarity with Quantum Mechanics, the 
variables with respect to which the brackets are calculated are relabelled as x, p instead of a, (3. 
Take the case where the bracket is Moyal; Then the work of Moyal showed that if ip{x, t) obeys the 
Schrodinger equation 

ih^ip = Hj>, (13) 

then the Wigner distribution function 

f°° - ipy 

f(x,p,t)= ip(x -y,t)ip k (x + y,t)e~dy (14) 



satisfies the quantum Liouville equation 





dt f(x,p,t) = {H(x,p),f(x,p,t)} MB , (15) 



where the bracket is the Moyal bracket. Baker [16| proved a significant converse theorem; that any 
real valued distribution function satisfying (|T5j ) together with a normalisation condition, could be ex- 
pressed in terms of a Wigner distribution function where the wave function satisfies the Schrodinger 
equation (|i~3"D 0. Baker's work suggests that one way to satisfy the equations (fl2|) is to take an ansatz 



of the form 

2iivpy 



Zk= / Xk(x-y,a)^/j k (x + y,a)e * dy (16) 



2 A more recent discussion of this view of quantum mechanics can be found in[p7|. 
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provided fa, k = 1 . . . 4 satisfy the Schrodinger equation 

H{x,—,a)fa{x,a) = ——ip k {x,a) (17) 

and Xk(x,a) satisfy the conjugate equation. This result follows from the observation that 

i 

-(H*z k -z k *H) 

= EEl-^'^T^^^X^)" / y m Xk(x-y,a)fa(x + y,a)e—dy 
-H <-> z k 

# (x + y,p + id x )xk{x - y, <r)^ fc (a; + y, a)e~^dy - H ^ z k 

-OQ 



2iirpy 

H(x + y, id x+y )Xk{x - y, a)^ k {x + y, a)e~^dy 



oo 

oo 



_ ; 7": / 



H(x-y,id x ^ y )xk(x-y,(r)fa(x + y,a)e A dy. (18) 

In a similar fashion the products Zj * z k can be calculated thanks to the following Lemmaf]. 

e py f(x)*e py 'g(x) = ry^lf-Aj^^/W^W 

■ u ml nl 

k 

= e^ + yy(x + y')g(x-y) (19) 

Then, by definition 

/OO rOO 2inpy 2iirpy' 

/ Xj(x-y,(r)fa(x + y,a)e-^-Xk(x-y',t)fa(x + y',a)e~^-dydy' (20) 
-oo J— oo 

This gives after integration over the variable y — y' 

f°° f / 2inp(y + y') 

z J *z k = N jk Xj {x-y-y ,a)fa{x + y + y ,a)e » d(y + y), (21) 



where Nj k is a (cr dependent) normalisation factor; 

POO 

Nj k = / Xk( x — m, a)ipj(x + u, a)du. (22) 



There remain the 6 = 3 + conjugate equations of constraint. One obvious way in which these can 
be satisfied is to choose Nj k = rijSj k . This then leaves one further equation to be solved, namely, 

d a H + {z u zi} + {z 2 , z 2 } + {z 3 , z 3 } + {z A , z 4 } = 0. (23) 

This equation determines the a dependence of H. However this procedure gives only a limited class 
of solutions, except in the 4-dimensional situation, which will be analysed again later. In order to 
solve the full 10-dimensional constraints in general, it is necessary to modify the ansatz to take into 
account the group theoretical structure of the self-duality relations in 8-dimensions. 



3 This was suggested by an unpublished result of Ian Strachan 
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6 Solution of 10-d equations 

Suppose 7 J , j = 1 ... 8 are the gamma matrices in 8-dimensions, which admit a real representation 
by 16 x 16 antisymmetric matrices [II| and ip is a 16 component spinor, which will be later subject 
to algebraic constraints. Take as ansatz 

POO _ 2iivpy 

X = i ip(x - y, cr)7 ip(x + y, aje^^dy, k = 1 . . . 8 

J — oo 

X 9 = i [ tp(x - y, a)^ 9 ip(x + y, a)e i ^dy, a = 9 (24) 



and leave x 9 as an unspecified anti-Hermitean matrix for the moment. By construction all X k , k = 
1 . . . 8, are real. There is nothing in this choice of which violates any symmetry principle; we are 
already dealing with a broken symmetry situation. Consider the six equations of constraint in (|T0|), 
which are identical to six of the seven equations of self-duality in 8-dimensions proposed in ||11|| . It 
is relatively easy to see that these last six will be satisfied automatically if if} is constrained to be 
of the form 

iP = P X P 2 P$ (25) 
where if} is arbitrary and Pi, P 2 , P3 are three mutually commuting projection operators given by 

Pi = J(1 + 7V7 5 7 7 )(1+7 2 7 4 7V) 

Pi = i(l + 7V7 8 7 5 )(l + 7 2 7 3 7 7 7 6 ) (26) 

p 3 = l (1 + 7 l 7 6 7V)(1 +727 5 7 3 78) _ 

This leaves the equations involving a dependence; the most conspicuous among which is 

d a X 9 = -({X 1 , X 2 } + {X 3 , X 4 } + {X 5 , X 6 } + {X 7 , X 8 }). (27) 
The other equations involving a derivatives and fl27D are then solved with the ansatz; 

7 9 = -(7 1 7 2 + 7 3 7 4 + 7 5 7 6 + 7 7 7 8 ), (28) 

if if} has an exponential a dependence. It might be thought that this method would lead to a general 
construction of a a independent solution, by choosing instead of (27) 



P 4 = 1(1 + 7 1 7 6 7 2 7 5 )(1 + 7 4 7 7 7 3 7 8 ), (29) 

and choosing as projection 

iP = P 1 P 2 P 4 i> (30) 

Then not only do the equations of constraint vanish, but so does the right hand side of ( p7[) and 
we should have, in principle, a solution to the 8-dimensional self dual equations. However, this is 
not so, as it may be verified that the product P1P2P4 = 0, i.e. the projection operator vanishes! By 
taking only two projection operators, P 2 and P3 in the construction it is found that only 5 of the 6 
equations of constraint vanish. The consistency of the system may then be restored by introducing 
another co-ordinate X°, which plays a similar role to X 9 . 
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7 Digression on the 4-d Moyal Nahm equations 

There is an alternative way to look at the 4-dimensional version of the first order equations, which 
is closer in spirit to the initial discussion in which the equations are regarded as an extension of 
Quantum Mechanics. This relationship is obscured in the 10-dimensional case on account of the 
equations of constraint. Consider the 4-dimensional Moyal-Nahm equations, where again X° is set 
to a constant 

d^X 1 + {X 2 , X 3 } = 
d a X 2 + {X 3 , X 1 } = 

d a X 3 + {X\ X 2 } = 0. (31) 

and asume that the functions ipj (x) belong to a real orthonormal basis over the real line. Then take 
as ansatz 

X 1 = J{E j -E i ){E i -E k )4 E i- E tt'' / ij)j(x - y)e^ip k {x + y)dy (32) 

J — oo 

together with cyclic replacements. The Ej are constants. Then on account of the result fl2l| ) the 
Moyal-Nahm equations are satisfied identically. This solution can be extended in an obvious manner 
by separating the functions into three classes according to the residue of the index mod 3; 

j=OC POO 

X° = J2 V^+i - E 3j )(E 3j - E 3j+2 )e^~ E ^ / tf; 3j+1 (x-y)e^^ 3j+2 (x + y)dy 
j=o J -°° 

j=0O roo 

X 1 = £ \J(E3 j+ 2 ~ E 3]+1 )(E 3j+1 - E 3j )e^~ Es > / ^ 3j+2 (x - y)e^^ 3j (x + y)dy 
j= o J -°° 

j = 0O /"OO 

X 2 = £ yj&j - E 3j+2 ){E 3j+2 - E 33+ M Ea] - E ' i3+l)u / ~ y)e^ip 3j+1 (x + y)dy 

j=Q J -°° 

The link with the above discussion on the analogy with the phase space formulation of Quantum 

Mechanics is as follows; When the correspondence p — > is implemented 

2n 

X 3 ^i(x)e- £lCT = ^(E, - E 3 )(E 3 - E 2 )e- E *° [°° ^(x - y)e yd ^ 2 (x + y)dy^ x {x) 

J —oo 

\E 1 - E 3 ){E 3 - E 2 )e- E *° / ^(ar - y)M^)M^ ~ v)dy 

J —oo 

{E x - E 3 )(E 3 - E 2 )M*>~ E2a (33) 



on account of the orthogonality relations. Similarly, exploiting the substitution y i— ► — y in the 
integral in (|32|) a similar equation for the adjoint appears; 



ij 2 (x)e E ^X 3 = -J(E 1 -E 3 )(E 3 -E 2 )e E ^ H V2 (a#i(x + y)e- yd ^ 2 (x - y)dy 

J —oo 

= - y J(E 1 -E 3 ){E 3 -E 2 )4> 1 (x)e El ° f°° ^ 2 {x - y)fa(x - y))dy 

J —oo 

= —\J(Ei — E 3 )(E 3 — E 2 )il)i(x)e Eia (34) 

This pair of equations can be combined into a Schrodinger like form with X3 playing the role of 
Hamiltonian. 
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7.1 Back to Matrix Models 

The treatment of the Moyal version of these first order Matrix Theory equations given here allows 
a neat return to the Matrix formulation in the following manner; define 

^ fc = EE Ak mn(t) r <Ux - y)M* + y)^dy (35) 

m n J oo 

where the coefficients A^ nn are a or r dependent elements of infinite-dimensional matrices labelled 
by the same index as X k . The use of the Lemma (|2l|) transforms the Moyal version of equations 
( |TUD into the infinite matrix version of the equations. This manipulation shows how closely the 
infinite matrix formulation is tied to the Moyal treatment. 



8 Conclusions 

This article has been written to examine the Moyal version of Matrix theory, in the belief that this 
formulation is the most appropriate for the discussion of the large N limit and for the investigation 
of parallels with Quantum Mechanics. The essence of the solution to the first order equations lies 
in the representation of the co-ordinates in terms of matrix elements of a spinor in 8-dimensions. 
The construction carries echoes of a similar representation of the 10-dimensional string in terms of 
bilinears |Tj| , but here a noteworthy feature of the construction is that it is intrinsically non-local. 



It would appear that instead of 'magic, mystery, or membrane' which have been suggested as a 
possible etymology the M in M-Theory really stands for Moyal! 
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